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O . Abstract 



We study the boundary free energy of the XXZ spin-1/2 chain subject to diagonal 
boundary fields. We first show that the representation for its finite Trotter number ap- 
proximant obtained by Gohmann, Bortz and Frahm is related to the partition function 
of the six-vertex model with reflecting ends. Building on the Tsuchiya determinant 
representation for the latter quantity we are able to take the infinite Trotter number 
iyy . limit. This yields a representation for the surface free energy which involves the solu- 

tion of the non-linear integral equation that governs the thermodynamics of the XXZ 
spin-1/2 chain subject to periodic boundary conditions. We show that this integral 
representation allows one to extract the low-T asymptotic behavior of the boundary 
magnetization at finite external magnetic field on the one hand and numerically plot 
this function on the other hand. 
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Among various models arizing in physics, one-dimensional models seem to play a very specific role. On the one 
hand, the reduced dimensionality of the system allows one to use various approximation techniques (correspon- 
dence with a Luttinger liquid [21] or a conformal field theory [9]) to be able to provide predictions for various 
quantities describing a given gapless quantum Hamiltonian in certain limiting cases: structure of the low-lying 
excitation, low-temperature behavior, long-distance asymptotic behavior of the correlation functions, etc. On the 
other hand, as it was originally observed by H. Bethe [5] on the example of the so-called XXX chain, there exists 
a rather large class of one-dimensional models that have the property of quantum integrability. In other words, one 
is able to characterize the eigenvectors and eigenvalues of the associated Hamiltonians with the help of solutions 
to certain algebraic equations, the so-called Bethe equations. The latter provide a very effective description of the 
spectrum of the model when the large volume limit L — > +oo is considered (ie when one deals with a model having 
a very large amount of pseudo-particles in its ground state). Such results allow one to identify the universality 
classes of various integrable models and, as such, check the predictions of the aforementioned approximation 
techniques, at least in some cases. For instance, using the Bethe Ansatz description of the spectrum, it was shown 
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to be possible to access to the 1/L corrections to the ground and excited states just above it, hence allowing one 
to identify the central charge and the various scaling dimensions [30, 31]. Such types of results gave explicit pre- 
dictions for the critical exponents governing the long-distance asymptotic behavior of two-point functions. The 
long-distance asymptotic behavior was then confirmed by extracting it directly from the exact representations for 
the correlators first at the free fermion points [38, 48, 52] and then for general interacting models [28, 36]. 

The question of the thermodynamics of Bethe Ansatz solvable models was first addressed by Yang & Yang 
[53]. These authors derived a non-linear integral equation whose solution allows one to compute the free energy 
of the so-called non-linear Schrodinger model [37] at finite temperature T and in the presence of an external 
chemical potential. The reasoning of Yang & Yang was raised to the level of a theorem by Dorlas, Lewis and 
Pule [13] with the help of large deviation techniques. The approach of Yang & Yang has been generalized to the 
study of the thermodynamics of the XXZ spin-1/2 chain simultaneously and independently by Gaudin [18] and 
Takahashi [45]. Their approach built on the so-called string conjecture [5] and allowed them to characterize the 
thermodynamics of the model in terms of a solution to an infinite hierarchy of non-linear integral equations. Since 
then, it has been applied with success to many other models and bears the name of the thermodynamic Bethe 
Ansatz. 

Another path to the study of thermodynamics of spin chains has been proposed by Koma [32, 33]. Building 
on the method for mapping quantum Hamiltonians in D-dimensions into models of classical statistical physics in 
D+l dimensions proposed by Suzuki [42], he argued that the computation of the partition function of the XXX and 
XXZ models in a magnetic field is equivalent to obtaining the largest eigenvalue of the transfer matrix associated 
with a specific inhomogeneous six-vertex model, the so-called quantum transfer matrix (QTM). Although Koma 
could not provide at the time a proper analytic framework for taking the so-called infinite Trotter number limit, he 
was able to carry out a numerical analysis along with an extrapolation to infinite Trotter numbers. Then Takahashi 
refined Koma's approach and was able to take the infinite Trotter number limit analytically. This led to a descrip- 
tion of the thermodynamics of the XYZ and XXZ spin- 1 /2 chains in terms of an infinite sequence of numbers 
that ought to be fixed numerically [46, 47]. Soon after, Kliimper [29] proposed an important simplification to the 
QTM-based approach. Namely, building on the method of non-linear integral equations [4], he proposed a way for 
sending the Trotter number N to infinity in the Bethe equations describing the largest eigenvalue of the QTM. He 
obtained the full description of the thermodynamics of the XYZ (and XXZ) spin- 1 /2 chains in terms of a single 
unknown function that satisfies a single non-linear integral equation. This function allows one to compute the 
free energy at any finite temperature T . Kliimper also proposed non-linear integral equations describing the sub- 
leading eigenvalues of the QTM what gave access to the correlation lengths. His results confirmed the conformal 
field theory-based predictions [1] for the low-T behavior of the free energy. 

All of the above results were obtained in the case of models subject to periodic boundary conditions. It so 
happens that the situation is definitely much less understood in the case of integrable models subject to other types 
of boundary conditions. Models such as the XXZ spin- 1 /2 chain subject to diagonal boundary fields 

M-i cosh (2 ) 

c^j + 2 +cosh tj (cr m o~ m+ \ + l)J+sinh tj coth ^ cTj+sinhT/cofh^-i- + _u/ \ (0.1) 

m=l 

have been solved through the coordinated Bethe Ansatz [3] and later by the algebraic Bethe Ansatz [41]. Above, 
cr x p , cr'p and cr z p are Pauli matrices acting on the Hilbert space V p =* C 2 attached to the p th site of the chain, n 
and are parameters characterizing the anisotropy and the boundary fields and J is an overall coupling constant 
fixing the energy scale. The analysis of the Bethe equations arizing in the diagonalization of the Hamiltonian (0.1) 
showed that the boundaries produce additional 0(1) contributions (ie ones that do not scale with the volume M) 
to the energies of the ground and excited states in respect to the values obtained for periodic boundary conditions 
[3]. 
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When turning the temperature on, one expects that changing the boundary conditions of the XXZ spin- 1/2 
chain will not alter the leading part of the model's free energy, ie the one scaling with the volume M. However, 
there will arize additional, volume-independent, contributions to the free energy. These constitute the so-called 
surface free energy. Various approximation techniques have been developed to characterize the surface free energy 
of the open XXZ spin- 1/2 chain (and also of more general, not necessarily exactly solvable models). Field theory 
based method led to predictions for the leading behavior at low-temperatures of the surface free energy in the 
vanishing external magnetic field limit [2, 7, 8, 14, 17]. 

The estimation of the surface free energy has also been considered in the framework of the thermodynamic 
Bethe Ansatz, what allowed one to make predictions in the full range of temperatures. Such calculations were first 
performed in the papers [11, 16, 54]. However, their results did not agree with the field-theory based results in the 
low-temperature limit. The origin of this discrepancy was elucidated in the papers [39, 50, 51] which re-considered 
the derivation of free energy of Bethe Ansatz solvable models, both in the periodic and the open boundary case. 
It was shown in [39, 51] that a careful derivation of the partition function leads to additional contributions, which 
are expressed in terms of Fredholm-determinants. These results were in agreement with the series of multiple 
integrals found previously in [12] in the case of integrable relativistic quantum field theories. The methods of 
[39, 5 1] have not yet been applied to the XXZ spin chain, but it is expected that analogous additional contributions 
are present in the case of the open spin chain as well, on top of the results given in [1 1, 16, 54]. 

However, all the aforementioned approaches to the calculation of the surface free energy where neither direct 
nor exact. A approach to obtain exact results for the surface free energy was initiated by Gohmann, Bortz and 
Frahm in [6]. The main result of that paper was a rigorous 1 representation for a finite Trotter number approximant 
of the boundary free energy. The latter was expressed as an expectation value of a large number of local operators 
forming the so-called finite temperature boundary operator. This expectation value was to be computed in respect 
to the eigenvector associated with the dominant eigenvalue of the quantum transfer matrix arizing in the description 
of the thermodynamics of the periodic XXZ spin- 1/2 chain. However, the representation given in [6] is rather 
implicit and it was not clear how to evaluate it at finite Trotter numbers, or how to take the infinite Trotter number 
limit. 

The present paper is devoted to overcoming these difficulties. Starting from Gohmann et al's result we show 
that one can, in fact, interpret the correlation function involving the finite temperature boundary operator as a 
specific case of the partition function of the six- vertex model with reflecting ends. Building on the techniques 
developed by Izergin [22] for the partition function of the six- vertex model with domain wall boundary conditions, 
Tsuchiya [49] showed that the former also admits a finite-size determinant representation. One cannot take the 
Trotter limit immediately on the level of Tsuchiya's determinant. We thus recast it in a form where the Trotter limit 
can be preformed easily. This yields the main result of this paper, namely an exact representation for the surface 
free energy of the spin- 1/2 XXZ chain subject to diagonal boundary fields. Using our exact representation, we 
obtain a simple integral representation for the boundary magnetization. This allows us to obtain the first terms of 
the low-temperature asymptotic behavior of the boundary magnetization. In particular, we show that we recover 
the results of [24, 26] for the T = case. We also use our representation to carry out numerical plots of this 
quantity in the massive regime of the model. 

This paper is organized as follows. In section 1, we briefly review the setting of the algebraic Bethe Ansatz 
framework for integrable models subject to the so-called diagonal boundary conditions. This allows us to set the 
quantum transfer matrix-based approach to the thermodynamics of the model and then review the representation 
for the finite Trotter number approximant of the surface free energy obtained by Gohmann, Bortz and Frahm. 
In section 2, building on a factorization of Gohmann et a/'s formula, we recast the aforementioned quantity in 
terms of Tsuchiya determinants. In section 3, we carry out several transformations on the formula obtained 

'modulo the strongly supported conjecture on the non-degeneracy of the quantum transfer matrix's largest eigenvalue and interchange- 
ability of the Trotter and infinite volume limits 
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in section 2, what ultimately allows us to take the infinite Trotter number limit. Finally, in section 4, we use 
our exact representation for the surface free energy so as to obtain an integral representation for the boundary 
magnetization at finite temperature. This integral representation allows us to extract the low-temperature behavior 
of the boundary magnetization in the gapless phase as well as to plot it, on the basis of numerical calculations, 
versus various parameters of the model such as the boundary magnetic field or the anisotropy. 



1 Preliminary definitions and the algebraic Bethe Ansatz framework 
1.1 Definition of the surface free energy 

One can diagonalize the Hamiltonian (0.1) within the algebraic Bethe Ansatz approach. For this purpose, one 
builds a one-parameter commutative family of operators, the so-called boundary transfer matrix [41]. The latter is 
expressed as a trace over a two-dimensional auxiliary space V a labelled by a roman index a of a product of bulk 
monodromy matrices T a (A), T a (A) and scalar solutions K*(A) to the so-called reflection equations [10]: 

r(A) = tx a [K + a {A)T a {A)K-{A)f a {A)] . (1.1) 

The bulk monodromy matrices are 2 x 2 matrices in the auxiliary space whose entries are operators acting on the 
quanta space b = ®^ =] V P of the chain. Thus, t(A) is also an operator on b = ®^ =] V p . The bulk monodrmoy 
matrices associated with the XXZ chain are built as ordered products of the 6- vertex type /^-matrix 2 



R(A) 



( sinh(/l + 77) 

sinh(/l) sinh(77) 

sinh(?7) smh(A) 

sinh(/l + 77) j 



(1.2) 



This /^-matrix can be seen as the matrix representation of an operator acting on a tensor product of two dimensional 
vector spaces C 2 ® C 2 . In the following, the subscripts a, b in R a b(A) mean that it acts non-trivially (ie not as the 
identity) solely on the tensor product V a ® Vb, where V a V\, ~ C 2 . Having introduced enough notations, we are 
now in position to write the explicit realisation for T a (A) and T a (A): 

T a (A) = R aM (A -&*).. .R*\ (A - ft) and % (A) = R la 01 + ft) . . . R Ma (A + ft,) . (1.3) 

There ft represent inhomogeneity parameters. We do stress that the roman index a refers to an auxiliary two- 
dimensional space whereas the indices 1 , . . . , M refer to the various quantum spaces V\ , . . . , Vm associated with 
the sites of the chain. As we have already specified, when these indices occur in a matrix, they label the spaces 
(auxiliary and quantum) where the latter acts non-trivially. 

Lastly, the expression for t(A) also involves the diagonal solutions of the reflection equations that have been 
first found by Cherednik [10] 

Kt (A) = K a (A + 77/2 ± ?7/2; ft) with K a {A; ft = | Sinh( ^J + ^ _ ^ J . (1.4) 

These /f-matrices can be checked to satisfy 

MJC(0)] =2sinh(ft)cosh(77) and tr B [^(0)] = 2sinh(ft) . (1.5) 



2 which, for further convenience, we chose to write in its polynomial normalization 
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sinh(w - v + rf) 
sinh(7/) 



-> + = 



sinh(w - v) 
sinh(?7) 



+ 



Figure 1 : The weights of the 6- vertex model. Here u is attached to the horizontal line and v to the vertical one. 

The /^-matrix given above satisfies a certain amount of properties. For instance, when the spectral parameter A is 
set to zero, they reduce to permutation operators 

R ab (0) = smh(r])P ab , (1.6) 

where P ab is the permutation operator on V a ® V b . They also fullfill a crossing relation 

o y l R t l 2 {A-rf)o y x =-R2 l {-X) . (1.7) 

Above, the superscript t\ in the /^-matrix refers to the transposition on the first space. The repeated use of the 
crossing relation results in the alternative representation for the boundary transfer matrix, 

t(A) = (-l) M tr a [K:(A)T a (A)K-(A)alT t a "(-A-rj)a y a } . (1.8) 

Further, bulding on the identities (1.6)-(1.5), one can show [10] that, in the homogeneous limit (& = 0, k = 
I,..., M), t(A) enjoys the properties 

, m ^a[K + a {0)}\x a [K-{0)} 2M /sinned 

r(0) = a - a - [sinh(77)] /M id and K = — -T(i) Ll=n . (1.9) 

There id stands for the identity operator on I) and "H is given by (0.1). As a consequence, 

[r(-P).r-'(0)f = e-?.(l + Or)) with fi = ^f™. (1.10) 

The above estimates allow one to write a Trotter limit-based representation for the partition function of the open 
XXZ spin- 1 /2 chain in a uniform external magnetic field: 



■h h N 

Z M = tri,..., M [e-£] with *Hh = <H - —Yt^k ■ (1-H) 



fc=i 

This representation takes the form 

M 

Z M = ^lirn o tr 1> ... (M |[T(- J 8/A0 ■ t~\0)] fje*^} . (1.12) 

a=l 

We stress that the subscript 1, . . . , M is there to indicate that the traces in (1. 12) are taken over all the quantum 
spaces V p , ie over the whole quantum space I) of the chain. Equation (1.12) admits a graphical representation 
given in Fig. 2. There, the summation over all possible ±-values of the intermediate lines in undercurrent and 
periodic boundary conditions are imposed in the vertical direction. The expression for the bulk (resp. boundary) 
weights is depicted in Fig. 1 (resp. Fig. 3). 
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-T] + p/N 

-fi/N 

-rj+fi/N 

-PIN 

-rj+P/N 

-PIN 



Figure 2: Graphical representation of the partition function (1.12) at zero overall magnetic field. The M vertical 
lines correspond to the quantum spaces of the spins, whereas the 2N horizontal lines joined by the curly end- 
curves correspond to the auxiliary space occurring in the definition (1.1) for the boundary transfer matrix. The 
bulk weights are represented on Fig. 1 whereas the boundary ones are represented on Fig. 3. Periodic boundary 
conditions are imposed in the vertical direction. 



+ sinh(-yS/A^ + £_) S~ sinh(/3/7V + £.) 



V2sinh<f_ V V2sinh<?_ 

+ ~~^\ _ sinh(-piN + 77 + £ + ) ~ _ sinhjpIN - ij + £+) 

V2 sinh £ + cosh(?7) . V2 sinh £ + cosh(?/) 

Figure 3: The boundary weights of the partition function depicted in Fig. 2. 
After some manipulations, as shown in [6], one is able to recast (1.12) as 



16=0 



(1.13) 



Note that t q (A) = tr[T^(A)] stands for the quantum transfer matrix associated with the quantum monodromy matrix 
arizing in the study of the XXZ spin- 1/2 periodic chain: 

^=^(^-^Witf-^.-.^«-^^tf-We* oJ = ( £g> ^ . (1.14) 

The quantum transfer matrix t q (A) acts on the tensor product of 2N two-dimensional spaces V ai ® • • • ® V U2N . The 
alternative expresssion (1.13) for the partition functions also involves the one-dimensional projector 

P ab {A) = Z?U)f* K~{A) . (1.15) 
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We remind that Tab is the permutation operator in V a <8> Vj,. 

It has been argued on the basis of thorough numerical computations, small ft analysis and free fermion point 
calculations [15, 19] that the eigenspace associated with its highest eigenvalue of the QTM for finite N is one 
dimensional. We shall build on this fact and denote by | *Po ) an associated eigenvector. We denote by Ao the 
so-called dominant eigenvalue associated to | ¥0 )■ Then, 

M ^Usinh 2/V (77)j \tr a [Ki (0)]tr Q [K- (0)]) ' <T |T > 

(1.16) 

Building on the rigorous approach to thermodynamic limits of partition function (see e.g. the monograph of Ruelle 
[40]), it has been argued in [33, 43] that, when computing the thermodynamic (M — > +00) limit of Zm, one can 
exchange the order of the N — » +00 limit with the infinite volume M — > +00 one. A consequence of such a fact 
would be that, when computing limM->+oo {(lnZ^)/M}, one can simply drop out all the terms in (1.16) that have 
been included in the . . . symbol. Indeed, for a finite Trotter number, all these terms will only produce corrections 
that are, for fixed N, exponentially small in M and thus vanish when M — > +00 and N is fixed. We shall not 
investigate this question further and simply assume that the exchangeability of limits does hold. 

The surface free energy is then defined as the limit of the difference between the free energy of the model 
subject to periodic and the one subject to open boundary conditions 

In virtue of the previous arguing, the surface free energy for the open spin- 1/2 XXZ chain is given by the below 
Trotter limit: 

.%t ,. r . _C (yo\Pa ia2 (-P/N)...P aiN _ iaiN (-/3/N)\y ) 

hm e t where e t = LJ 1N — (1.18) 



T 

W->+oo 



<T | ¥o> • {tr a [K (0)] tr a [K~ (0)] /2) 



The representation for e"^ as the expectation value (1.18) constitutes the main result obtained in reference [6]. 
There, it was also observed that (1.15) is indeed a one-dimensional projector. This can be easily seen as soon as 
one observes that the transpose of the permutation matrix is as a one-dimensional projector 

Kl = [I + >«l +)b + \~ )a\ ~ )b\ ■ [( + |«< + \b + < - l«< - \bl (1.19) 

where | ± ) a is the canonical spin up/down basis in V a . This structure allows one to factorize the representation for 
the surface free energy at finite Trotter number. Indeed, setting 

|V> - (l+>a,l + >* + |->ai|->fl2)®-"®(l + WMl + W + |-W-il-W) ( L2 °) 
(v| - « + | ai < + la 2 + (-L,(-L)® •••®(( + k v _ 1 ( + k v +<-L 2w _,(-U,v) (1.21) 



_lsurf 

one recasts e t as 



. a _ ( 2 \ N < l^i i-PIN) . . . K + UlN l {-PIN) I v > 



v tr [tf + (0)]tro[*-(0)]/ <^o|^o> 

x ( v \K~ { {-pIN) . . . K~ 2N l {-PIN) I To > (1.22) 

The representation ( 1 .22) will constitute the starting point of our analysis. Indeed, we will establish that each of 
the two expectation values occurring in the numerator in (1.22) are related to the partition function of the six- vertex 
model with reflecting ends. However, we first need to discuss in more details the construction of the eigenvectors 
of the quantum transfer matrix. 
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1.2 Eigenvectors of the QTM 



The QTM can be diagonalized by means of the algebraic Bethe Ansatz [32, 33]. Indeed, the eigenstates of 
f = u& [T l H^) ] are build as a repetitive action of B operators on the pseudo-vacuum 



\~)a 



|0> = |+> B1 ®|->a 2 ® 

More precisely, one has that 

t\t)-B(A L )...B(A l )\0) = [sinh(77)] 27V T | {A k \^) ■ B (A L ) . . . B (A^^) , 

where 



(1.23) 



(1.24) 



k=\ 



sinh(£ - ^) 

L 

k=i 



sinh (rj) 



sinh(£ - ^ + rj) t sinh(£ + /3/iV - 77) sinh(£ - /3/AO 



sinh(£ - A k ) 



sinh (77) 



N 



(1.25) 



is the eigenvalue of t q (^) [sinh(?7)] 2N associated with this choice of A^s. The parameters A k are subject to the 
Bethe Ansatz equations 



-1 - e~T 



n 

k=\ 



sinh(/l p - Ak + rj) 
sinh(/L -A k -rj) 



smh(A p + BIN - r/) siah(A p - BIN) 



sinh(i p - BIN + 77) sinh(i p + B/N) 



(1.26) 



Numerical investigations, analysis at the free fermion point and for B small indicate that the dominant eigenvalue 
is given by the choice L - N, this for any value of h [20]. What changes, however, is the actual distribution of the 
roots. These are located on the purely imaginary axis when h = and occupy regions with a more complicated 
shape as soon as h + 0. 

From now on, we focus on the solution {A^ describing the dominant eigenvalue. We shall also agree that 
I ¥0 ) refers to the Bethe Ansatz-issued eigenvector | ¥0 ) = FI^li B (^a) 1 ). As proposed in [29], it is useful to 
introduce the function closely related to the exponent of the counting function 



_h 

e T 



1 — r sinh(o» - Ak + rj) 
11 sinh(w -A k -rj) 



sinh(w + B/N - rj) sinh(w - B/N) 



N 



(1.27) 



sinh(<x> - B/N + rj) sinh(ct» + B/N) 

The "a function is OT-periodic, bounded when %(oj) — > +00 and such that it has, in the case of generic parameters, 

• an N th -order pole at co = -B/N , 

• an Af th -order pole at a> = B/N - rj , 

• N simple poles at oj = Ak + rj, k - 1 , . . . , N. 

It thus follows that 1 +"0 (a>) has 3N zeroes. N of these are, by construction, the Bethe roots, ie 

1 + a(A k ) = , k = l,...,N. (1.28) 

Numerical analysis and calculations at the free fermion point indicate that the roots for the ground state of the 
QTM can all be encircled by a unique loop that is moreover ,/V independent and such that any additional root 
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to 1 + ~a (a>) = is located outside of this loop [20]. It appears, on a numerical analysis basis, that all these other 
roots accumulate around ±n. In the following, we shall build our analysis on the same assumption. We do stress 
that all of the above considerations should be accommodated so as to be consistent with the natural in periodicity 
of the functions involved. Putting all these information together allows one to conclude that the function"? solves 
the non-linear integral equation 



_ h 
In a((o) = -— + Nln 



sinh(w + B/N + rj) sinh(a> - (3 IN) 
sinh(w - J/N + tj) sinh(w + J/N) 



£e^(jU)-n)}n[l + aQi)]-^. (1.29) 



where we have set 9(A) = /In [sinh(?7 - A)/ sinh(?7 + A)]. The contour encircles the Bethe roots A\,. . . , An, the 
pole at -B/N but not any other singularity of the integrand. We also remind that, in the derivation of the non-linear 
integral equation, one has to make use of the fact that ln[l +"a(o>)] has a zero monodromy around (the contour 
encloses the N th order pole at -B/N and the N simple zeroes at A\, . . . ,A^ of 1 + a(a>)). 



1.3 The norm of Bethe states 

Under the normalization of the /^-matrix that we have chosen, the "norm" of the eigenstates of the QTM admits 
the determinant representation [34] : 

N 

<0|]~[C(^)-]~[B(^)|0> = m ^Jl.[ 5 (Aj,p/N)5(Aj,/3/N-T])] N dety[/+Z] . 

/=i i= l J= l 1 a( }) ' n sinh(i a - A b ) 

a+b 

Here, det>y [/ + K] is the Fredholm determinant of the trace class integral operator I + K acting on L 2 Ctf) with the 
integral kernel 



K{co,o)') = - ^ ^ where K(A) = G'(A) 



i sinh(2?7) 



(1 + a{o)'))2n ' sinh(/l - rf) sinh(^ + 77) 

Finally, the "norm" formula also involves the shorthand notation 
s(A,fi) = sinh(/l - ji) sinh(/l + jj.) . 



(1.30) 



(1.31) 



2 Rewriting the expectation values 

In this section we first recast the expectation values occurring in (1.22) in terms of a semi-homogeneous limit 
of the partition function of the six- vertex model with reflecting ends. As observed by Tsuchiya [49], the latter 
admits a determinant representation. A direct calculation of the semi-homogeneous limit through the L'Hopital 
rule is however unadapted for our purpose. Hence, in the second part of this section, we apply the Cauchy 
determinant factorization [23, 28] so as to recast the Tsuchiya determinant into a form allowing us to take the 
semi-homogeneous limit in an elegant manner. This ultimately leads us to obtain a relatively simple explicit 
representation for the finite Trotter number approximant of the surface free energy. 

2.1 Relation with the partition function 
Proposition 2.1 The expectation values 

T~ = (v\K- ai (-p/N)...K- 2N l {-/3/N)B(A 1 )...B(A N )\0) (2.1) 
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and 

T + = < |C Ui) . . . C 01*) < (-/3/AO . . . K 2N J-J3/N) | v > (2.2) 
occurring in (1.22) can be recast as 

<f ~ = e~w(Q | C~(-B/N) . . .C'(-B/N) | 0) and T + = e^(0 | ® + (-B/N) . . . S + (-B/N) | 0) . (2.3) 

N terms iV terms 

There C~, resp. S + , refer to entries of the boundary monodromy matrices of -, resp. +, type: 

<U~ a (A) = T a (A) K~(A) T a (A) = ( ^ (2.4) 

=Tt(A)[K;(A)y°ThA)= ™) , (2.5) 

T a (w) =/? a jv(«-^iv)...^ai(M-^i) and T fl («) = /?i fl (w + /li) . . .Rno(u + A N ) . (2.6) 
Lastly, we agree upon 

|0)=|~)i®---®l~);v and |0)=|+)i®-'-®l+)y. (2-7) 
with | ±)# being the canonical spin up/down basis in the space Vj — C 2 . 

Proof — 

Using that B (A) = k {+ \ T*(A) \ -) k , we get 

r~ = ( v | ® (0 1 (-/}/#) . . . ^ {-BIN) R^ - /3/W) Ui - )S/AO . . . 

. ..^(-Ai -B/WRuMi -jS/AOe*^ . ..^ 2 y-% -B/N)R Nai {A N -jB/JV)e&°*|0>® |0) 
= < v | ® (0 1 <^ - /3/W) . . . - jS/JV) ^ (-)8/A0 /W. 0*1 - jS/JV) . . . 

N 



/W> (An ~ BIN) . . . K~ (-BIN) R^M - BIN) . . . R Na] (A N - B/N) J~[ (e*"*) • 1 > ® 1 0) (2.8) 



By using (2.6), we reconstruct monodromy matrices in which the Bethe roots {A^ for the ground state of the 
quantum transfer matrix play the role of inhomogeneities and -B/N is interpreted as the spectral parameter. This 
leads to 

r~ = e"$<v|® (0| (-yS/A^)^^^ (->S/A^) • • • ^ (-^/A^)^;, (-/5/A^) • |0>® |0) 

It now remains to compute the partial scalar products involving the even spaces a^, k = 1, . . . , N. This can be 
done thanks to the identity 

U-t<-U<-| + a 2 a + M + l) -< a 2 f(A)\-) a2k = a2 a-\Ta 2k -M) (2-9) 
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Therefore, 




( < - 1 ...«,<- 1 ® (0 1 {-PIN) . . . {-PIN) | + ), 



...| + > fll ®|0) 



- e _ 5r (0 1 C~(-p/N) . . . C~(-P/N) 1 0) . (2. 10) 



Very similar steps can be applied so as to re-write T + . Indeed, one can express C {A) = *(- 1 T^{A) \ +\ what 



T + = (0| ® (0| ^(-^i -PimRu^Ui -PIN) . ..^Hi -plN)R lai {A\ -PIN)c^ . . . 

■ ■ ■ R i N (- A N ~ PIN) R Ncn Un ~ PIN) e*^< (-0/N) . . . K^_ t {-p/N) | v > ® 1 0) 
= e^<0|<g> (0| ^(-Ai -/3/AO . ..S^Ni-Ax-fi/mRlaurMl -PI N ) ■ ■ -Rnow-Mn ~ PIN) 

('PIN) ...RiaMi-PIN)... R Nai Un - PIN) < {-PIN) | v > ® 1 0) 
= e$< 1 ® (0| (-j8/iV) 7^ (-y3/A0 (->8/A0 Taw-i (~P/N) 



In the case of ;F + , one takes the partial scalar products in respect to the odd spaces a2k-i, k - 1, . . . , N. This is 
done by means of the identity 



J + \T a2k _, (AW:^ {A)(\- ) m _, | - ) aik + | + ) a2k _ t | + ) a2k ) = [K 2k {A)p k 7% {A) | + ) cl2k ■ (2. 12) 



T + = e® ai J - 1 . . . a ( - 1 ® (0| [u- 1N {-piN)f w . . [u- 2 {-p/N)f 2 | + ) a2N . . . | + > fl2 |0) 

- (0 | S + {-plN) . . . ® + {-p/N) | 0) . (2. 13) 



The scalar product in expression (2.1) is depicted on Fig. 4, whereas the first expression in (2.3) is shown in 
Fig. 5. Note that we have used the crossing relation so as to recast part of the weights into a canonical form. Note 
also that these two figures are related through a reflection across the North-East diagonal. The proof of proposition 

2.1 corresponds to an algebraic verification of this symmetry. 

2.2 Alternative factorization of the Tsuchiya Determinant 

It follows from the previous observations that the two factors J r± can be identified with the partition function of 
the six-vertex model with reflecting ends. The latter has been computed, in terms of a determinant, by Tsuchiya 
[49]. More precisely, let C be given by (2.4), then the aforementioned partition function is defined by 



leads to 



. . . T a a 2 2 {-PIN) T ai {-PIN) K+ {-PIN) | v > ® 1 0) (2.11) 



Therefore, 



Thus proving the second identity. 




(2.14) 



By using the relations provided in [26], it is readily seen that 



(0|£ + (^)...S + (£v)|0) = f] 



a=l 



sinh(2£, + 2?]) 
sinh(2£,) 



) 



(2.15) 



n 



A\ 



A N 



-> + 



-> + 



-> + 



+ 



+ + + + + + + + 

J3/N PIN p/N PIN 

J]- PIN rj-p/N rj-p/N tj-P/N 



Figure 4: Partition function of the 6- vertex model with a reflecting end. This graph corresponds to the expression 
(2. 1): the horizontal lines represent the action of the B-operators on the reference state, whereas the open boundary 
conditions at the top correspond to the action of the local ^-matrices. 

The Tsuchiya determinant representation for Z,n ({ftdf ', Wfc}^ £-) reads 



J! j sinh(£_ + A a ) sinh(2^)| • n {^a, h) • s(ft + 77, A h )) 
Zh({£^-MU-) = - z — ^[/V] 



N 

n 

a<b 



(2.16) 



fl {s(/lfo, A a ) sinhfe - ft) sinh(ft, + ft + 77)} 



s(/l,/i) has been defined in (1.31) whereas the entries of the matrix are given by 
_ sinh(?7) 

The antisymmetry properties of the determinant along with its invariance under the transformations 



{A a 



,N 



{o- a A a } 



and 



{ftJi ^ - + l)77/2}f for any e a , cr a e {±} 



(2.17) 



(2.18) 



ensures that the singularities present in the denominator of the right hand side effectively cancel out. Taking the 
homogeneous limit £ a — » ~P/N recasts (2.16) in terms of a Wronskian. Such a representation is however not 
adapted for our goals. Notwithstanding, by using the Cauchy determinant factorization techniques [23, 28], one 
can provide a representation for Z,n that cancels out explicitly the apparent singularities at ^ = ft and A a = A},. 



Lemma 2.1 The determinant of the matrix N can be decomposed into the below product of determinants 



det w [N] = f] 



a=l 



1 



sinh(2/l a ) 



• det w [ 



1 



sinh(ft - Aj) sinh(ft + Aj + rj) 



— ] • det w [6jk + U jk ] 



(2.19) 
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+ + + + 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



-r] + pIN 



-PIN 



-r] + pIN 



-PIN 



-r] + pIN 



-PIN 



-r] + pIN 



-PIN 



-An 



Figure 5: Partition function of the 6- vertex model with a reflecting end. This graph corresponds to the first 
expression in (2.3): The horizontal lines together with the boundary weight correspond to the action of the C 
operators, and the fixed boundary conditions on the top and the bottom describe the two reference states (0| and 
10). 



Where the entries of the matrix Uj k read 
sinh(2/l<; + rj) 



sir±(A k + Aj) sinh(i ; - -A k + r]) *r . 

Ft sinhOl* - A a ) 

0=1 



J! SWh(A k -Aa-Tj) N 
<;=! Y\ " + k 



0=1 



sinh(i fl + A k + n) 



N 

n 

o=l 



sinhfe - A k ) sinhfe + A k + n) 
sinhfe + A k ) sinh(£, - A k + rf) 



(2.20) 



Proof — 

Given the Cauchy matrix 



1 



sinh(& - Aj) sinh(& + Aj + 77) 



(2.21) 
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Its inverse has entries 

]1 [ sinh(ft - Aj) smh(Aj + ft + 77) smh(ft - A a ) sinh(i a + ft + 77)) 

[C '] 1 



(ft ^) * + ^ + V) { sinh(ft - ft) sinh(ft + ft + 77)} ft j sinh(^ - /IJ sinh(/l; + ^ + 77)) 
It is readily seen that 



fl=l 0=1 



with 



N 

N FI sinh(ft - Aj) smh(Aj + ft + 77) 

X^Vv^ 5 * < 2 - 22) 

P =i fi sinh(/l; - sinh(/b + A a + 77) 



2V 

/v . , t , FI sinh(ft, - A a ) sinh(i a + £ p + 77) 

c = V sinh ^ £=± (2 23) 

>* Zisinh(^-^)sinh(^ + ^ + 77)5(^,^)5(^ + 77,^) n . , v 

p - ! FI smh(ft - ft.) smh(ft. + ft + 77) 

£1=1 



Now, note that, on the one hand, the /7r-periodicity of the below integrand leads to 



N 



(£— sinh(7/) sinh(2oj + 77) pr sinh(oj - A a ) sinh(/l a + + 77) (? ^ ^ 



J 4in sinh(a> - Aj) sinh(<x> + Aj + T])5(a>, A k )s(oj + 77, A k ) 1 _1 sinh(ft - sinh(<jj + ft + 77) 

OT-strip 



On the other hand, the integral can be taken by computing the residues located in the OT-periodic strip. These two 
observations show that 

_ 1 sinh(/l ft - /W sinh(A Q + 4 + 77) 

]k ~ Jk smh(2A k ) fl* 1 sinh(ft - A k ) sinh(A k + ft + 77) 

+ 1 pr sinh(/lfc + /l fl ) sinh(^ - A a -if) ^ 

sinh(/l/t + /I,) sinh(/l, - ^ + 77) sinhCl/l^) A A sinh(ft ( + Ak) sinh(ft - A k + 77) 

CZ=1 

Hence, after replacing 5^ in (2.22) by (2.25), pulling out the pre-factors and carrying out a similarity transforma- 
tion, we obtain the desired representation. ■ 

Proposition 2.2 Let {A a }^ be a solution of the Bethe equations. Then, the homogeneous limit (ft — > —B/N) of the 
partition function admits the below representation: 

z4-BiNtM a )U.) ^^tM^^lM ) . fi smh( ^;;^7 ) 

v ' l a= \ I smh(2A a ) I 1J sinh(/l fl + /l 6 ) 

iv 

x P] [ sinh(/l fl - B/N) sinh(?7 - ^ - B/N)] ■ det N [6 jk + U jk ] . (2.26) 

a=l 
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where the matrix Uj k reads 

N 

_h . II smh(A k - Aa + rj) N 
jj k= ~ e T sinn ( 2 ^ +JV a=i pi sinh(/t a + A k ) ^ 

Jk sinh(A k + Aj) smh(Aj -A k + rj) *t . \\ sinh(^ a + A k + n) ' 

EI sinh(Ai - A a ) a ~ l 

0=1 

Proof — 

Using that the parameters {A a }^ satisfy the Bethe equations, it is readily seen that U j k = Uj k . The rest follows 
after straightforward algebra. ■ 

As follows from proposition 2.2, most of the apparent singularities of Z,n can be canceled out by means of 
the Cauchy determinant factorization. Yet, in (2.26) there are still apparent singularities at A k = -Aj. These are of 
course compensated by the zeroes of the determinant det/v [6j k + U j k ]. However, their presence is problematic in 
respect to taking the infinite Trotter number limit. We now factor out these zeroes explicitly. 

Lemma 2.2 The below factorisation holds 

det N [S jk + U jk ] = f[[l + oM,)]* • f^fif • • ( 2 - 2 8) 

p=i ' ' 

where^a has been defined in (1.27) and 

-in-k 



+ °° C +C ° f/(^2/t+l)l 2k+1 d 2k+ 

r N {{A a] ») = 2 j) 2 2 , + 1 J • 11 ^ 7i"y 



+oo 



- X 9 Z aT , n^ l) — • ,129) 



ylm) 



Above we agree upon to n +i = to\. Also, the integrands contain the function 

N 

sinh(/l a + to - rj) sinh(/l a - to - rj) 



j (to) = e T \ \ (2.30) 

i | sinh(/l a + w + n) sinh(/l a - to + rj) 



as well as the kernel 



-~ -e / smn(zw + rj) i-r sinn^ + w ; sinntw - /t a + r]) 

U(aj,u> ) = . (2.31) 

sinh(a»' + to) sinh(<x> - to' + rj) i_| sinh(&/ - A a ) sinh(/t a + to' + rj) 

Finally, for any p,^ D • • • D ^ are encased contours such that ^ k , for k = 1 , . . . , p enlaces the roots A\ A^ 

but not the ones that are shifted by ±rj. 

We would like to stress that the series (2.29) might not be convergent. However, the function eT^iUa)^) is 
well-defined. More precisely, one can take 

N . _ . I 



f][l + a(-A p )]-i • (|-^|) 4 • det w [S jk + U jk ] , (2.32) 
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as its definition. It is readily seen by inspection of the zeroes and poles of the right hand side of this equation, 
that the left hand side does not have singularities at Aj = -A^. More details about the precise mechanism for this 
definition are given in the proof below. 

Proof — 

Let \k\ be small enough, so that the logarithm of the determinant admits the series expansion 



+°° K 2n+l N +°° N 

In detjv [6jk + Uj k ] = J] ^— j-y J] U PlP2 . . . U P2n+m - ^ — ^ u pm ■ ■ ■ u P 2n Pi ( 2 - 33 ) 

n=0 Pl,-,P2n+\ n=\ Pl,-,P2n 

= 1 =1 

We first derive a recurrence equation that allows one to compute the various traces. Namely, given any pi,p2,P3 € 
H 1 ; N J and agreeing upon the shorthand notation A a b - A a + A\, and A a b - A a - A\, one gets 

N 

n 11 sinh(^ + 77) N - 

U P1P2 U P2P3 = -e"T sinh(2^ 3 + 77)— jj ■ . ,7 : S pm ■ < 2 - 34 ) 



£i nsinh(i P3fl ) l a=\^Ua P3 +v) 



p 3 c 

a=l 
*P3 



where, 



N 

n . mi , > 11 sinh(^ P2a + 77) N , 

V - sinh(2^p, + 77) a=i 12 r-r smh(A ap2 ) 



£=i sinh(i P2/ , 3 ) sinh(^ P2P3 + 77) sinh(^ lP2 ) anh(yl WK +77) ^ sinh(/} ^ sinh(/t ap2 + 77) 

a=l 
*Pl 

C du> - sinh(2(jj + 77) si 

j 2ot sinh(o> + A P3 )smh(a> - A P3 + 77) sinh(tu + A Pl ) sinh(/L, - a> + 77) A J si 



sinh(/l a + (J) sinh(a> - A a + 77) 

■ sinh(o> + /1 P3 ) sinh(tu - /l P3 + 77) sinh(<jj + A pi ) sinh(/t pi - <x> + 77) A J sinh(<x> - A a ) sinh(/l a + a> + 77) 

/v 

n sinh(i pia ) Af -j , 
_ s 1 ggi p| sinhQ^p, - 77) 

PljP3 sinh(2^ Pl + 77) n _ 1 1 &mh(A p . a - ?]) ' l ' ) 
P[ n sinh(i pifl ) «=i pifl 

a=l 

Above, the contour ^ encircles all solutions {A a }^, but not the ones that are sifted by ±77, ie the sets [A a ± 77}^. 
Thus, 



2_j U PIP2 U P2P3 ~ /</W<W 3 + <I> sinh(6> + 



e rry^^aOsinh^/l^ +77) dw 
/lp 3 ) sinh(<x> - /l P3 + 77) 2in 



N 

J! sinh(i P3a + 77) a, • , ,7 v 

Ff - — (2-36) 

nsinh01 P3fl ) a=i sinh (^3 +? 7) 

a=l 
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When dealing with the contour integral-based term, one can compute the remaining sums over A pi , A P3 , A P4 . . . , A Pn 
by a similar contour integral provided that one successively uses an encased contour ^D-o^,. This choice 
of contour ensures that the poles at to p = -a) p+ \ do not contribute to the value of the integral. Ultimately, one gets 

N r n — d n u> N 

V U PuP2 ...U Pn , Pl = (p r\u(co p ,co p+ i)- — ^- + V f(A Pl )U PuP2 ...U Pn _ 2 , Pl . (2.37) 

U~,Pn ^_ J _^P=l K ' PU-,Pn-2 



Where we agree upon to n+ i - co\. 

The induction can be solved and leads, in the even case, to 

N n-l „ fc 2( "~*L 2(H ~ k) d0J N — 

2 U Pl , P 2---U P2n , Pl =J] <b [f(a>2(.n-k))] Y[u(a> p ,co p+l )-Y\ + £ [f(A Pl )f . (2.38) 

Pl,---,P2n k=Q ^ p=\ p=\ pi = \ 

■■0^2(n-*) 

Whereas, in the odd case, 

2 U PUP2 ■ ■ ■ U P2n+uPl = £ (p [/( W2(n _, )+1 )] ] ] U(co p , to p+l ) • HI it + Tj \f&pJYU PlPl ■ 

Pi,-jn n+ \ k=o ^ P =i P =i P i=i 



(2.39) 



Yet, 

Pi=l ^ a=l 

It then remains to observe that 

f(0) = [a(0)] 2 and /U pi ) = -a(-^,) . (2.41) 
what leads to 

N n n , 2(n-t)+l j2(n-*:)+l j 



2 U Pl,P2 ■■ -Upturn =Yt 9 [/( W 2(n^) + l)] f] U&P'Up+l)- . ,2(»-fc) + l + 2'["° (0) ] 

■■■,P2 n+ l fe=0 p=l 



2n+l 



■■0^2(n-«+l 



(2.42) 



Finally after inserting the two formulae in the trace-like expansion for the determinant, weighting by appropriate 
powers of k and using 



ln(l + x) - ln(l - x) = V ^^x" = 2 V , (2.43) 



n=l n=0 

we are led to the representation 

N 



det, [s jk + ^] = ]> + ■ (r^ij ■ ^ Wl) ■ (2 - 44) 
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where 



"Z"* j> Zh^ n^M^I. (2-45) 

Since the kernel U and the function / are uniformly bounded on the sequence of integration contours, it is 
readily seen that the above series convegres for |/c| small enough. One thus has an equality between holomorphic 
functions of k, in some neighborhood of k = 

e f«({-Uf) = + a:-A p )]-5 • ( |~^ ) 4 • det^v [fy + rffyt] . (2.46) 

p=i ' ' 

The function in the rhs is an analytic function on the closed unit disk. This thus ensures that the function k h-> 
e^N ) is analytic on the closed unit disk. It can thus be continued analytically from a neighborhood of k = 
-where it is definitely licit to define it through the formula (2.45)- up to k = 1. If the series (2.45) is convergent at 
k — 1, then we get the representation (2.28) with Tn ) given by (2.29). Else, the representation (2.28) still 

holds but the definition of exp [T~n ) } is to be understood in the sense of analytic continuation from an open 

neighborhood of k = up to k = 1 of the function k i-» exp yj^ ({/l a }^) } with ({/l a }^) defined, for \k\ small 
enough, by the series (2.45). ■ 



2.3 The representation at finite Trotter number 

By putting together all of the previously obtained formulae, one arrives to the below representation for the finite 
Trotter number approximant of the surface free energy: 

fjurt m t~t sinh(£_ + /l a )sinh(f + + A a ) t^I sinh(?y) sinh(-2/3/A0 
= 6 T H sinh(£_)sinh(£ + ) 1 = ) sinh(2i a + rj) ' \j sinh(2i fl ) 



~f sinh(i fl +A b + T]) f Uatb sinh(i a - A b )\ [ \\ N a=x t 1 + \ ( 1 + a(0) * 5 



n^inh^-^+T,) I ft^^^ J ln^sinh(^ + ^)J\l-a(0) y 

\ sinh(2?7) j dety[/ + ]fl ' 



n a'(^)M^) 
a=i 



3 Taking the infinite Trotter number limit 

3.1 Rewriting of the double products 

The representation (2.47) constitutes a good starting point for taking the infinite Trotter number limit. For this, as 
it is customary in the QTM approach, one should represent all simple and double products over the Bethe roots 
for the largest eigenvalue of the QTM in terms of contour integrals involving the function "a defined in (1.27). 
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Once such a representation is obtained, the infinite Trotter number limit can be easily taken. The purpose of the 
proposition below is to provide such a contour integral representation. 

Proposition 3.1 The below products admit the alternative representations 

f[*bn.-*) = f\lM ■ i^M^dmif . n (— (jlj^mh)" . i «-* f 

U [J "Mil 1 sinhto-W/V) / yisinhfo + ^+AW)/ V(l +e -f) 2 ' 

exp|-2A^^^ln[l + a"(ai)]coth(6j+/3/iV) - ^^-coth'(w'-w)ln[l+'a(a;)]ln[l+'a(oj')]| ■ 

■if V tf'cV 

(3.1) 

A/so 

at 2 N I +"a(-/l ) f Td 

[1 sinh(i a + A b ) = [sinh(-2/3/A0f ' • f~[ exp - 2W (f)-^ ln[l coth (to - PIN) 



x exp j (j)^ (|)^-cothV + w')ln[l + a(w)]ln[l + a(a/)]j ■ (3.2) 

A sinh(2^) 1 + a(0) f r dto\ 

a-l ^ 

A sinh(^ + A fe + ??) = | sinhfo-^/AQ )* 2 A f sinhfo + A a + /3/AQ ) N 
}J l smh(A a -A b + r ] ) \ sinhfa) J H \ sinh(77 - ^ - jS/W) J 

x exp | - 2N ln (1 +"a(a>)) [coth (a> - P/N + if) - coth (w + p/N + 77)] | 

( f da> dco' _ . . _ , . / sinh(a» + a>' + 77) \ 1 

xexp (f)— — hi (l + a(w) In (1 + 0(0') RA/ln — -J ^ . (3.4) 

[ J 2?7r 2j7t \ sinh(a» - a)' + rf)j ) 

The contour ^ is the one occurring in the non-linear integral equation (1.29). It is such that the poles of the 
integrand at to = ±p/N - r/, resp. to' + to + rj - 0, are all located outside of 'it?, or'&X'tf, depending on the integral 
of interest. 

We do stress that the encased contour c €' D corresponds to a small deformation of ^ such that the poles in 
to' - ±co are all always outside of 16". 

Proof — 
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The representation for the fourth product is the easiest to obtain. Namely, 



a,b=\ 



smh(a> + A a + rf) \ a'(a>) dco sr-\ / smh(-fi/N + A a + rf) 



Z A l n s i nh (^a + Ab + if) _y i smin,g> t- A a + rp \ u yoj) uuj ^ yi i 
11 sinh(/l fl -A b + rf) ~ "I sinh(a> - + 77)/ ' 1 + a" (w) 2ot + 4^ ° I sinh(-/3/7V - /l a + 77) y 



j sinh(a> + <j/ + 77) \ 'a'(cu) a'(a>') dwdo/ /\ /sinh(w -J3/N + rj)\ a'(a>) dco 



^ ' \sinhu-j - a/ + rpj 1 -'a(oj) 1 + a(a/) (2/>) 2 " T "' \ sinh(w + (3/N + ;/) / 1 +a(w)2/7r 

y l ( smh(+fi/N + A a + 7 1 ) \ 2 fa / sinh(77 - 2/3/AQ \ f / smhfrj + 77 ~/3/AQ \ dcu 
^ n \ sinh(-j8/W -A a + rf)J \ sinh(?7) / J \ sinh(w + fi/N + rf) ) ' 1 + a"(<y) 2ot 



(3.5) 



Formula (3.4) then follows after an integration by parts. Doing so is licit in as much as In' [1 + !x(co)] has a 
vanishing monodromy along ^ . Note that the contour ^€ is chosen precisely so that the poles of the integrand in 
respect to co (resp. co') at co ± co' + 77 = mod[/7r] are all located outside of c €, this for any co' e <?f (resp. co e c €). 
Next, we consider 

N 

S\ («) = ^ In [sinh(/i a + At + «)] . (3.6) 

a,6=l 

It is then easy to see that, for u small enough, 

d<y (cu) v-i a \~A a - u) 



S\(u) = ) j (b—coth(A a + oj + u)- 1/ + > NcotHA a -]3 N + u) - ) J —^ > , 

jziJ 2 ™ l + a(<u) £f ^l + a(-^-M) 

fda> fd<j/ , "a'((t») 'a'(co') „ IT f da> , , ^, 1TN "a ' (a>) 

= (Dtt- CD 77- coth(w + a/ + if) • 1/ \' + 2W (h — coth(w + K-/3/A0 • - ' - 
J 2in J 2m l + o (a>) 1 + a(a» ) J 2j^ 1 + a (a>) 

, a'(B/N-u) ^ 'a'(-^-w) 

+ Ar 2 coth( M -2/3/A^) - N ^- > — . (3.7) 

l + a(fi/N-u) j^l + a(-A a -u) 

Above, the contour c €' is such that it does not encircle the points -co - u, with 00 € c to. Taking the integral over u 
and carrying out integrations by parts, we obtain 

C* 'dec) C* ddc^' dfx^ 
= m— (D — coth'(aj+w'+M)-ln[l+a(w)]ln[l+a(a>')] - 2AT (T) — coth(w + « -/3/A^>ln[l+a(w)] 
J 2in J 2in J Lin 

N 

+ Ci + A^ 2 ln[sinh(w-2/3/A^)] + ATln[l + a(J3/N - u)] + V ln[l +"a(-A a - m)] . (3.8) 

a=i 

There Ci is an integration constant that ought to be fixed. Taking u — > +00 in the original representation (3.6), we 
get that <Si(w) - A^ 2 [w-ln(2)] + 2N 2Za=i + o(l). The same limit can be taken on the level of the representation 
(3.8). Indeed, one can always send u — > +00 by deforming it (and the contours ^ , c €' if necessary) from u - 
along a path that keeps the properties of the contours entering in the double integral unaltered. Then, it remains to 
use that, for any fixed z,"a(z - u) - e~T(l + o(l)), whereas 

f _ do* £ 'a'(co) dco 
-(()ln[l + a(a>)]-— = (D oj- — )=--—■ ■ — - = fi + YA a . (3.9) 
J 2in J 1 + a(co) 2in ^ 
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These, in the u — > +00 limit, lead to 

N 



Si(u) = N 2 (u -In 2) + 2N^ A a + 2Nln[\ + e~r] + C\ + o(l) . (3.10) 



Thus, in order that the alternative representation (3.8) produces the correct large u asymptotics, one has to set 
C\ - -2AHn[l + e~r]. Further, setting u - in (3.8) and using that"o(/3/A0 = leads to equation (3.2). 
We now establish the factorization formula (3.1). For this, we introduce 

•S2 (u) = > In [sinh(/l a - Ab + u)] so that sinh(i a - Ab) = lim < — \ . (3.11) 

A »->o Uinh^fw)-' 

Then, 

Q6J (o>) V a + u > 



S' 2 (u) = yd)^cotHA a -u + u).-^- + y^Vcoth(i fl+ ^ + M ) + Y-^-^ 

^ J 2ot 1 + a ( w ) z_i *-f 1 + a(A„ + 

fda; fda/ , 0' (gj) aV) y ?(^+«) 
- CD CD cothi^j - <i) + m) • ~ ~ + > - 



(Dt- ([) — coth(a/ - 6> + 11) • 1/ L, ' + > — ^ + Af 2 coth(a) 

J 2m J 2in 1 + a (a>) 1 + a(a>') ^ 1 + a (A a + u) 

+ N £ coth (a; + u +J3/N) - coth (co + fi/N - u) ] ■ * ~? . - N *^f , *~ U \ . (3.12) 
J 2in y J l + a(cu) \ + a(-fi/N-u) 

Here, the contour c €' is such that to - u lies outside of c €' for any u e Thus, taking the anti-derivative in respect 
to m and then carrying out integrations by parts, one gets 



S 2 (u) = - <£^- (£^-coth'(w'-w+M)-ln[l+a(a>')]ln[l+o(w)] + ^ln[l+aU a + u)] + N 2 ln[sinh(w)] 



N 



j2in J 2in 



+ C 2 + N]n[l + a(-/3/N-u)] - N £[ coth (a) + u +J3/N) + coth (co +/3/N - u)] -In [1 • 

(3.13) 

The integration constant C2 can be fixed by carrying out much the same reasoning as for Si(u). One gets that 
C2 = -2AHn[l + e"r]. It then remains to take the exponent and compute the u — > limit as in (3.11). This can be 
done by observing that 



N 



iim , n l + a(^ M ) i j-j 
k->0 I A 1 sinh(w) J 1 1 

a=l a=l 



as well as 



M r — >> t^t sinh(-/L + 77 — /3/AO f sinh(-n) sinh(w) l w r n w 

sinh^fw) 1 + a(-u-BN)] ~ e~r 5: — 2 — 7 y - L_ y -L • sinh(-28/A01 

1 v H ;] u->o 11 sinh(-i Q -?7-/3/A^) lsinh(77 - 2^/AO sinh(-w)/ L H 1 

N 



N * A sinh(,-^-^) ( sinh(,) J* f f 15 
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N 

Finally, we compute the last product. Setting S3 (u) = 2 In [sinh(2/i a + w)] we get that : 

0=1 

S' 3 {u) = £ coth (2co + u) tt 'ff ^ - 1 • - ° i~7 /2 ^ + # coth (-2/W + u) . (3. 16) 

J \ + a(co)2in 2 1 + a(-u/2) 

Hence, after an integration in respect to u, 

S 3 (w) = -2^coth(2w + M )ln[l+T(w)]^ + In [1 +~a(-u/2) ] + AHnsinh [(-20 /N + u)] + C 3 . (3.17) 

<€ 

The integration constant is fixed by comparing the u — > +00 asymptotics of the two representations for S$(u). One 
gets, on the one hand, 

N 

S 3 (w) = 2Y,l a + N(u-\a2) + o(l) , (3.18) 

0=1 

whereas, on the other hand 

N 

S 3 (u) = 2(Y^ A a + P) + ln(l+e"7) + N(u - In 2 - 2/3/N) + C 3 + o(l) . (3.19) 

a=l 

This implies that C3 = — ln(l + e~r). ■ 

3.2 A smooth representation 

Inserting the previous formulae into (2.47) leads to 



e t 



I / sinh^sinh^^/AQ f «p{-[g(g + ) + gtf-)]/r} 1 + ^ , ( , ?1 

= ~^~T7 0o/AA ■ r /0 : j = = f =T ex P \ 2Tn V A *h ) + J • ( 3 - 2 °) 

\smh(?7-2/3/A^)smh(2^)/ J[ _ rj (0 )i2 det»|/ + JC| 1 V ' ' 



Above, we have set 



Vl - RO)] 2 de% [/ + A 1 ] 

= —N In ( sMl ^ " ;f /A ° ) - Ttfg In (1 + a(-£» + r (fin (1 + a(w)) coth + (3.21) 
v sinh(f) 7 j 27/r 



where we agree that 

_. f 1 if - £ is located inside of the contour 

op — i „ , . • (3.22) 

s [0 otherwise 

Furthermore, we have introduced 

1 = 2 j^)^- In (1 +'a{co)) [ coth (2w) + coth (2co + 77) ] 

- 2N (j)^- In (1 + a~ (a>)) [coth (w + /3/AO - coth (a» - /3/AO + coth - /3/N + 77) - coth (w + PIN + 77)] 



- coth'(w-w') + d^d,,/ In — - coth'(w+w ) 

L sinh(w - to' + 77) J 



+ CpTT- (D ^-ln(l +"a(w))ln(l + a(w')) 

J 2j/T J 277T 



(3.23) 
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and we insist that the contour arizing in the second line is such that the points in/2, -rj/2 mod[/7r] are outside 
of <g. 

Finally, the function Tn (Ma)f ) is given as in (2.29) with the sole difference that now, the kernel U, is given 

by 



U(a>, oJ) 



-e t sinh(2a/ + 77) 



sinh(a» + a>') sinh(<x> - 10' - rf) 



sinh(a>' - /3/N) sinh(w' + /3/N + 77) 



sinh(a>' + /3/N) sinh(w' -/3/N + 77) 



exp I - ^jr- In [1 + (r) ][ coth(r + to') + coth(a/ - r) - coth(r + to' + 77) - coth(<x/ - t + 77)] | , 



(3.24) 



whereas the function / (to) is recast in the form 



sinh(w - /3/N - 77) sinh(w + (3/N + 77) 



sinh(w + /3/N - 77) smh(a> - fi/N + rj) 



exp I - ^ > ~2 r ' m [1 + a (t) ][coth(r + to - 77) + coth(w - r - 77) - coth(r + w + 77) - coth(a> - r + 77)j| . 

(3.25) 

The contour ^/ is such that given any to', to e where ^ refers to any of the encasted contours introduced in 
lemma 2.2, the points 

±oj', ±(a)' + rj) ±(a/-Tj) (3.26) 
are not surrounded by c tojj. The latter loop encircles however the points {A a }^ as well as the origin. 

3.3 Representation for the surface free energy 

In order to take the infinite Trotter number limit, one should first send N — > +00 on the level of the non-linear 
integral equation for the function"?. There N only appears in the driving term which has a well defined N — > +00 
limit. It thus appears highly plausible that"? — > a where a is the solution to 



h 2 J sinh(?7) 
hxa(.a>) = -- + 



{ coth(w + 77) - coth(w)} + ^(co -fi) • ln[l + a(/j)] • ^ . (3.27) 



Once the question of the limit of the function a is settled, it is not a problem to send N — > +00 in the above 
formulae. One gets 

/ surf = + S(£_)-2/3r[coth(77) - coth(277)] - 2T-T - T I 

T 1 
+ -ln[l - <4 g (0)] + rin(deV[/ + Z]) - rin[l + e"?] (3.28) 



In this representation, we agree upon 

S(£) = r/3coth(£) - T5 r in [1 + a(-£)] + T (^) In [1 + a(w)] 



da» 

coth(w + £)• — . 

Z77T 



(3.29) 
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with 8^ as defined in (3.22). a reg (0) corresponds to the N — > +00 limit of ^T(0). We do stress that the Trotter limit 
and the co — > limits do not commute for'a(w). In fact, a reg (0) corresponds to the regular part of the function a(o>) 
at co = 0, viz. 

lna reg (0) = limflnaM + — sinh(77)coth(w)l = -- + — coshfo) + (f) tffr) -In [1 + afci)] • = . (3.30) 
Furthermore, one has 

I = £ In (1 + a (a))) [2 coth (2w) + 2 coth (2cj + 77) - 4/3 coth' (co) + 4/3 coth' (co + 77)] ^ 

r sinh(a) + + 77) n 

flcAi* ln — ) - coHi' (co-co') - coth' (co+co') 

L smh(o) - co' + 77) J 



jTdcu f dco' , , ,„ , ,^T„ „ , r sinh(a) + co' + 77) - 



+ (Dtt- (p — ln (1 + a (a/)) ln (1 + a(a/)) 

J 27/T J 277T 



(3.31) 



Also, the kernel K is defined as in (1.30) with the sole difference that the function a ought to be replaced by 
the function a defined by (3.27). Finally, T is given by 



™ r ^[f(co 2k+1 )f k M d 2* + l w 



* =( Md..o% +1 "= k - :1 



[/("*)] 



n-k 



2 * j2£ 



d co 

x 4> z ,„ • n^^)^' (3 - 32) 



The kernel £/ defining the function ;F is given by 

-e _ 7 sinh(2a»' + 77) , , , 

co') = — —A- £ • exp - 2/3 coth(w') - coth(a/ + 77) 

sinh(a» + co') sinh(w - co' - if) 

x exp j - (|) jr- In [1 + a (t) ] • [ coth(r + co') + coth(o/ - t) - coth(r + co' + 77) - coth(a»' - r + 77)] j , 

(3.33) 

and the function / (o») reads 

f (co) = e T exp I - 2/3 [coth(w - 77) - coth(w + 77)] J 

exp j - ^jr- In [1 + o(r) ] • [ coth(r + co' - 77) + coth(a>' -t-tj) - coth(r + co' + 77) - coth(a/ -r + ?7)]| . 

(3.34) 

The encased contours ^ are as defined in lemma 2.2 whereas the contour ^7 is such that given any a/, w € 'tfp, 
p € N, the points ±o/, ±(o/ + 77) ± (a/ - 77) are not encircled by ^7. The latter however encircles the region where 
the numbers A\,...,A.n condensate and, in particular, the origin. 
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4 The boundary magnetization 



4.1 An integral representation at finite temperature 

It follows from general considerations that the boundary magnetization can be obtained from partial derivatives 
of the partition function. More precisely, one has 

sinh 2 (£_) d r 

(<ti)t-m = ^ • -rr lnZ M . (4.1) 

Hence, using that 

^M-expj- — T + 0(-)J (4.2) 

and assuming that it is licit to exchange the M — > oo limit with differentiation, we get that, in the thermodynamic 
limit 

sinh 2 (^) d .... 

^1 fT =~ T . 7 , \ ■ T^/surf • (4.3) 

J sinh(?7) ot, - 

In fact, solely the function S(£~) defined in (3.29) gives a non-zero contribution to (4.3): 

rsinh 2 (g-) ggjg-g-) ^ sinh 2 (^) X Ml + oH 

\0"i)r = 1 + • + r • — . (4.4) 

7sinh(77) l + ct(-£_) J sinh(?/) J sinh 2 (a> + f_) 2ot 

4.2 Numerics 

Building on the previous integral representation, we present numerical calculation issued plots for the boundary 
magnetization. We do not attempt a full exploration of the parameter space; our aim is to demonstrate that our 
results can be easily implemented by suitable computer programs and that the known T — > limit is reproduced 
numerically. Therefore we constrain our investigations to one of the possible domains in parameter space, namely 
the so-called massive regime (cosh(?7) > 1, rj € R) of the chain 3 . 

We choose to express quantities in units of / what amounts to setting / = 1 in the numerical analysis. 
Also, we found it convenient to recast the parametrization of the boundary magnetic fields in (0.1) as -hi = 
2sinh?7Cofh£ ± . 

We remind that, in the massive regime, it is more practical to consider the rotated counting function a m (A) = 
a(iA). The latter satisfies 

h 2 sinh 2 77 f sinh2?7log(l + a m (u))) dto 

log a m (A) + + (7) • — . (-45-) 

T T sm(A) sm(A - irj) J sm(A — a> + irj) sin(/l - oj - irj) 2n \ ■ > 

There, the integration contour ^ consists of two intervals: 

^ = [-ia - n/2, -ia + n/2] U [ia - n/2, ia + n/2] , 

since the right and left lateral contours cancel out in virtue of the ^--periodicity of the integrands. Finally, the 
parameter a < rj/2 is a real number which has to be chosen large enough so that all the roots parametrizing the 



3 The reason for this name is that the model has a finite gap between the (possibly doubly degenerate) ground state and the first excited 
states at sufficiently small external magnetic field. The excitations become however massless for sufficiently strong values of h 
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dominant eigenvalue of the QTM lie inside the contour. We have solved the integral equation (4.5) by the iteration 
method. 

Within these new notations, one gets the representation for the boundary spin (p\) (4.4) 

/ z\ _ i _ T sinh2 £- X" log(l + a,»(a>)) du _ sinh 2 ^_ ia' m {i£-) 

sinhrj f sinV-^-) ' ^ * " ' sinh 17 1 + a m (^_) ( 46 > 

We remind that 5^ = 1 if lies inside the contour and = otherwise. Since it is not convenient to evaluate 
derivatives on the numerical level, we have rather used the expression below for a numeric evaluation of the last 
term in (4.6) 

ia'ag-) sinh 2 77 sinh(-2£_ + 77) f 4 sinh 2?] sinh(2(£_ + ioj)) dco 

8 t - (J) log(l + a m (a>)) • — (a n\ 

a m (j£_) T [cosh(-2£_ + 77) - cosh 77] 2 J 6 [cosh(2(£_ + ico)) - cosh(2?7)] 2 2n^ n 

We numerically evaluated the boundary magnetization as a function of the magnetic fields for different values 
of the temperature and for 77 = 1.5. The choice a = 0.9rj/2 appeared appropriate for all the values of T and h 
that we have considered. In Fig. 6 we consider one specific point in the parameter space (77 = 1.5, T = 1, h = 2, 
h7 = 0) and show that indeed there is an interval a* < a < rj/2 where the numerical result for the boundary 
magnetization does not depend on a. 

Our plots are shown in Fig. 7, in what concerns (o~j) as a function of the bulk magnetic field at h~ b = 0, and 
Fig. 8 in what concerns (cr^) as a function of the boundary magnetic field at h = 0. We have also added the plots 
corresponding to the zero-temperature limit. The latter have been extracted by using the T = representation 
obtained in [24, 25, 26]. 

There arize two values of the external magnetic field in the description of the ground state of the Hamiltonian 
(0.1). At ht = h7 = and for h > Ii^ = 4(1 + cosh 77) the ground state is completely polarized and (cr?) = 1. 
Then for h c ) < h < h c ), with 

h ( £ = 4 sinh 77 V (4.8) 
cosh(«?7) 



the model is massless. Then at ht - h, = 0, the boundary magnetization admits the integral representation: 

(a\) = 1 - f cot(v-777/2)g(v)dv (4.9) 



-9 

There the function h solves the linear integral equation 



,^ _ f sinh(277) ■ g(r) dr _ d_ , sinh(?7) > 

J sin(A - t + irj) sin(/l - r - irj) 2n ds Wsin(/l + s + irj/2) sin(/l + s - ij]/2)>\s=o ' 
-1 

Note that the endpoint of integration q e [0;n/2] is defined as the solution to the equation 60(17) - 0, where 
the dressed energy eo is given by the solution to the integral equation: 

f sinh(277) dr 27sinh 2 (77) 

e ^ ~ ~^~71 , • \ ■ ,1 r-eo(r)-— = h . . — — — . (4.11) 

J sm(/l - t + 777) sm(A -t — itj) 2n sm(A + itj/2) sm(A - 777/2) 



26 



0.37 
0.36 
0.35 
0.34 
0.33 
0.32 
0.31 
0.3 
0.29 
0.28 
0.27 
0.26 



N p = 20 
N p = 100 



0.45 0.5 0.55 0.6 0.65 0.7 0.75 

a 



Figure 6: In this figure we demonstrate the dependence of the numerical results on the parameter a, which de- 
scribes the position of the integration contours. Here we consider one example with 77 = 1,5, T - 1 and h — 2. 
The integral has been discretized into N p = 20 and N p - 100 points. The theoretical requirement is that a < rj/2 
and a should be big enough so that the contour encircles all the Bethe roots of the Quantum Transfer Matrix. 
In the present case the first singularities appear around a « 0.47. Note that numerical errors arise also around 
a = 77/2 = 0.75, however their magnitude decreases with growing N p . The numerical value of the boundary 
magnetization is (cr^) = 0.356912, and it changes less than 10~ 6 in the middle regime, even with N p - 20. 

One has q = for h = f$ and q = jt/2 for h = h„. At h = h„, the model starts to become massive and, for any 
< h < hjj, one should set q = n/2 in (4.10) and (4.11). In such a case, the integral equation (4.10) becomes 
explicitly solvable via Fourier transformation. The boundary magnetization is then given by (4.9) with q = n/2. 
In this case, one can compute it in a closed form, this for any value of the boundary magnetic field h~ [24] : 

<°1> = 1 + 2 d - ^ E ( ;_^ r)2 where -± = — • sinh^ . (4.12) 

In Fig. 9 we plot results for the general situation with both the bulk and the boundary magnetic field present 
(here we fix the temperature to T = 1). 

In Fig. 10 we also plotted the boundary magnetization as a function of the anisotropy parameter rj, in the 
case of the ground state (T = 0) and also in a finite temperature situation (T = 0.5). In the latter case our data 
is limited to smaller values of q, because our simple iteration algorithm loses its good convergence properties at 
higher values of the ratio cosh(?7)/T. 



4.3 Low-r expansion of the boundary dependent part 

Proposition 4.1 In the massless regime of the chain q = —i£ with £ € ] ; n [, the boundary field dependent 
integral given by (3.29) admits the low-T asymptotic behavior 

8(g) = r/3coth(a + J cothO* + £ + ;<r/2)e(A) • ^ 

C f 

nT 2 sinh(2,f + it) A 

l2s Q (q) smn(q + g + ig/2) sinh(g - g - i£/2) 

where 

Q = [-q;q] U r(-i^/2 if < -5(g) < £/2 
= [-q;q] otherwise 
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Figure 7: The boundary magnetization as a function of the bulk field h, for different fixed temperatures (77 = 1.5 
and h7 = 0). The zero-temperature result is given by equation (4.9). The properties of the ground state change at 
the two critical values hl r - 2.6585 and h 2 rr - 13.410. 




Figure 8: The boundary magnetization as a function of the boundary field h b , for different fixed temperatures 
(rj = 1.5 and h = 0). The zero-temperature result is given by equation (4.12). 
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(a) Boundary magnetization in the ground state (T = 0) 



(b) Finite temperature results (T = 0.5) 



Figure 10: The boundary magnetization as a function of the anisotropy parameter rj, for different boundary mag- 
netic fields. The bulk field is set to h = 0. The zero-temperature result depicted in Fig. 10(a) is calculated using 
equation (4. 12). 

Furthermore, we denote by T(z) a small counterclockwise loop around z. 



Proof — 

In order to obtain the low-J expansion of the boundary field dependent part of the surface free energy, we first 
need to obtain the first few terms of the low-J asymptotic expansion of the function a. In the massless regime, 
within the parametrization r\ = -ig with < £ < n, the non-linear integral equation satisfied by this function takes 
the form 



ei)(co - i£/2) £ du r 



(4.15) 



Here we agree upon 
e Q (A) - h - 



27sin 2 (0 



sinh(i + i£/2) sinh(/l - i£/2) 



and 



K(ju) 



sin(20 



sinh(ju + ig) sinh(ju - i£) 



(4.16) 



The integration contour ^ can be chosen as depicted in Fig. 11. There a < min(^/2, (n-£)/2), and is large enough 
so that the contour 'rf encircles all the roots describing the dominant eigenvalue. 

We assume that we are in the anti-ferromagnetic regime of the chain meaning that the magnetic field is not 
sufficiently strong enough so as to polarize the chain to the ferromagnetic state, namely 



27 sin 2 (0 7 
0<h<h c with h c = — = — — = 8/cos 2 (<f/2) 

sin 2 (<T/2) 



(4.17) 



In such a situation, the function eo(A) admits two (symmetric) zeroes +qo on R. Moreover it satisfies 

e °\]- qo -. qo [ <0 and e m\[- qo m] >0 - < 4 - 18 ) 

We would like to analyse the behavior of the solution a(co) at low temperatures. There, it appears that the dominant 
contribution at T — » (up to 0(T°°) corrections) will stem from the line R + i£/2. It is apparent from the previous 
discussion that the contour ^ can be readily deformed towards this line for £ € ] ; n/2 ]. 



29 



% = -R + ia 



% = R - ia 



Figure 1 1 : Contour ^ for the non-linear integral equation for a in the massless regime. 

Indeed, suppose that < £ < n/2. One expects that, at low-temperatures, a(to) = 0(T°°) whenever a> belongs 
to the region encircled by the curve ^£ that moreover lies strictly in the lower-half plane. The function of a> defined 
by the integral 



tfV-^Ml + aC//)]^ (4.19) 



can be analytically continued from co lying inside of ^ up to a> e R + i£/2. One can even deform the integration 
contour ^ up to R + i£/2. Then, agreeing upon 

e(A) = -T In [a(A + if/2)] for leR, (4.20) 

we are lead to the below integral equation 



s(A) 



= e Q (A) + T J ^-K(A-fi)-ln[l +e _£ ^] - T J K(A - n + i£/2 + ia) ■ In [1 + aQi - ia)] . (4.21) 



One should, in fact, choose a in such a way that a(pi - ia) = 0{T°°) uniformly in n e R, this in the (L 1 n L°°)(R) 
sense. This condition should be checked a posteriori once that a has been computed. If its holds, this means 
that the low-J analysis of the non-linear integral equation reduces to the one of an equation quite close to the 
Yang- Yang equation 

K(co - n) ■ In [1 + e"— ] • ^ + 0(T°°) . (4.22) 

R 

Although it is unclear to us how such a deformation procedure would work for f € ] 7r/2 ; n [, we shall nonethe- 
less work under the hypothesis that this can be done. In other words, we shall assume that equation (4.22) holds 
for the whole regime < f < n and that the power-law part of the T — > asymptotics contribution stemming from 
integrations of analytic functions versus In [ 1 + a(A)\ along ^ always issues from a deformation of the contour ^ 
to R + i^/2. 

The equation (4.22) resembles to the Yang-Yang equation [53]. Techniques for extracting the low-J asymp- 
totic behavior of its solution are well known [15, 29]. A rigorous approach to the extraction of the \ow-T asymp- 
totic expansion of the solution to the Yang- Yang equation has been given in [35]. The setting of this last paper 
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can be directly applied here when n/2 < £ < n, as then K{A) < for A e R. For < £ < 7r/2, the change in 
the sign of Lieb's kernel K breaks certain properties used in [35]. However, one can still compute the asymptotic 
expansion on a formal level, by assuming the a priori existence thereof. Independently of the value of £ we shall 
give a formal derivation of the asymptotic expansion. Hence, we assume that, just as e$, e admits two roots ~q^ 
on R. Likewise, we shall take for granted that 

e l] ?<-);?<+)[< whereas s |R\[?h ;?<+>] > • (4.23) 

It follows from (4.21), by slightly shifting a if necessary, that e(A) is holomorphic in some open neighborhood of 
the real axis. The size of this neighborhood does not depend on T. Thus, according to the results of lemma A.l, 
one has 



In other words, the function s solves the integral equation 



?<-) 

There, the 0(r 4 ) is in (L 1 n L°°)(R). We assume that e admits the low-J asymptotic expansion 

e(A) = £Q (A) + Tei(A) + T 2 s 2 (A) + 0(r 3 ) (4.26) 

and that (f (±) — » ±^r in the T — > limit. It is easy to see that, under such assumptions, £o solves the linear integral 
equation 

i 

so (A) + | K(A - u) so (u) — = eo(X) with q fixed by the condition Eo(±q) - . (4.27) 
J 2/r 

In other words, as one could have expected, eo is to be identified with the dressed energy of the particle/hole type 
excitations above the ground state of an open XXZ spin- 1/2 chain at finite magnetic field h. Note that, when 
n/2 < £ < n, one can use the techniques developed in [35] so as to prove the unique solvability of (4.27) for eo 
and q. Equations (4.26) and (4.25) allow one to fix the dependence of the endpoints <f (±) on T. Namely, 

= s^) = s (±q) + (§<*> + q)s' (±q) + T £l (±q) + 0(T 2 + (q {±) + qf) . (4.28) 

Hence (q^ + q) goes to zero at least as T (this if £\{±q) + 0, otherwise it goes to zero even faster, ie at least 
as T 2 ). Inserting the above expansion (4.28) into the linear integral equation (4.25) and keeping terms that are at 
most a 0(r 2 ), we get 

?<+) q 

f K{co-n)s{n)^- = f K{^-n)e{^ + ^ -q)K{u-^) E -^ - (q<~> +^ (£U _5<->)f<0 + 0((q^hq) 2 ) 
J 2n J 2n 2n 2n 

?<-) -q 

(4.29) 
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The first two terms vanish since s(g^) = 0. Hence, the corrections stemming from the displacement of the 
endpoints in respect to ±q produce at most 0(T 2 ) corrections in (4.24). As a consequence s\ - 0. This means 
that, in fact, the corrections stemming from the displacement of the endpoints are at least a 0(T 2 ). Hence, as 
follows from (4.29), they generate 0(T 4 ) corrections in (4.24). After rearranging all of the expansions, we get that 

s(A) = s Q {A) + -^-\R{A,q) + R(A, -q)\ + 0(T 4 ) , (4.30) 

where R(A,/j) is the kernel of the resolvent operator to / + K/(2n), ie (I - R/(2n)) ■ (I + K/(2n)) = I. 

In order to access to the low-J asymptotic behavior of S(£) given in (3.29), just as in our analysis of the 
non-linear integral equation for the function a, we first deform the upper part of the contour to R + if /2 and, if 
necessary, slightly shift upwards the lower part of the contour. Upon applying lemma A. 1 , this leads to 



e(-;j/2-g C dA 

S(f) = r/?coth(f) - T8An[\+e t ] + coth(A + £ + i£/2)e(A)— 

J 2in 



;r 2 r 2 f coth(^+f + /f/2) _ coth(g<-> + f + if/2) \ 4 

2ot6 1 s>(qM) £'(?(")) J 1 J ' 1 ^ 



One can simplify the logarithm using that 

%[e(-i£/2 - f)] < for < -3(f) < f/2 
!R|>(-if /2 - £)] > otherwise 



(4.32) 



Also, the endpoints q l ~ ±) can be replaced by ±q since the former deviated from the latter by 0(T 2 ) corrections and 
one also has that £(§^) - 0. This leads to 

£(f) = Tfico±(£) + ri A (f) £( ~^ 2 ~ a + J cothOl + f + if/2) s{A)^- 

-i 

t:T 2 



+ i- 



\2e' Q {q) 



{ coth(<7 + f + if /2) + co\h{q - f - if /2)} + 0(T 4 ) . (4.33) 



Where 1^ represents the indicator function of the set A = {z € C : < -3(z) < f/2}. It is clear on the level of 
the expansion (4.33) that the second term can be recast as a contour integral leading to (4.13). ■ 

We can now apply the previous proposition so as to obtain the low-temperature expansion for the boundary 
magnetization. Indeed, inserting the low-J asymptotic expansion of S in the reconstruction formula (4.3) for (cr^) 
one gets 

sinh 2 (f_) C s (A) dA 



smh^ r 
x l/ /isin(f) J 



sinh 2 (i + f_ + if /2) 2ot 

, 7i sinh 2 (4L) d sinh(2£_ + if) , 

+ T 2 — ^ — + 0(r 4 ) (4 34) 

12/sin(£)£o(<7) <9f- sinh(# + f_ + if /2) sinh(# - f_ - if /2) 

One can check that the T = limit in this above expansion does reproduce the expression obtained in [26] 
at finite magnetic field h. Indeed, let us recast the first line of (4.34) in the language of that paper. One can first 



32 



carry out an integration by parts. The boundary terms coming from an integration along [-#;#] vanish since 
£o(±q) = whereas the contour integral r(-if/2 - f_), should it be present, has no boundaries. Thus 



r sinh(f )g (i) ^ _ _ C tq ^ sinh ^_) . coth ( A + ^_ + . dA (4.35) 

J sinh 2 (i + f _ + it IT\ J 



sinh^i + f_ + if /2) 
One then has that 

sinh(f_) cosh(^ + f_ + if/2) = cosh(f_) sinh(i + f_ + if/2) - sinh(/l + if /2) . (4.36) 

Clearly, the first term in the r.h.s. of (4.36) will only have a vanishing contribution to the integral so that 

sinh(f_) r e' (A) sinh(/l + if/2) 
) J 

It is easy to check that e' (A) solves the linear integral equation 



sinh(f_) f sinhQl + if /2) 
= 1 - J2^T) J sinh(J + f- + if/2) • ^ " (437) 



e&M) + f ^(i-ju)enO")— = -27r/sin(f)— ( ). (4.38) 

° v J v w oW 2tt vt <9AVsinh(A + if/2)siiih(^-if/2)' v ' 

The so-called inhomogeneous density of Bethe roots in the ground state p(/l;f) which was used in [26], satisfies 
the linear integral equation 



g 

p(A;f) + f K(A-fi)p(A;^ = — r^^- ... - ... ■ (4.39) 
J 2n n sinh(/l - f ) sinh cosh(?/)(/{ - f - if ) 



Thence 



1^/2 = , , (4.40) 
2;r7sin(f) 



and, upon substitution, we get 

inh(^ u- f if/2) 0t L __ : 



/ \ , f sinh(g-)sinhQt + if/2) a ,,,,,, 

<cri> r=0 = 1 - — ttt-; — — p(A;f;) dA . (4.41) 

J sir 



which is precisely the representation found in [26]. 



Conclusion 

In this paper we have studied the so-called boundary magnetization of an XXZ spin- 1/2 chain subject to diagonal 
boundary fields on each of its boundaries. Starting from a representation for its finite Trotter number approximant 
obtained by Gohmann, Bortz and Frahm [6] we have recast it as a product of partition functions of the six-vertex 
model with reflecting ends. Using the determinant representations of the latter partition functions obtained by 
Tsuchiya [49] along with the Cauchy determinant factorization trick [23, 28] we have been able to recast the 
resulting expression into a form that allowed us to take the infinite Trotter number limit explicitly. We then 
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applied our result to the computation of the boundary magnetization at finite temperature. In such a way, we were 
able to check that the zero temperature limit of our result does reproduce the known integral representations for 
(cTj } and also to draw some curves for the latter quantity in the massive regime. 

A natural continuation of our study would be to address the question of building an effective approach to the 
computation of the correlation functions in models subject to diagonal boundary conditions at finite temperature. 
Such a work would provide one with the boundary analog of the representations obtained in [20] and an extension 
of the representations for the correlation functions at T = OK obtained in [26, 27]. 

It would also be interesting to confront the predictions issuing from the thermodynamic Bethe Ansatz with 
our results. The full boundary free energy has not yet been obtained in the latter approach. More precisely, 
certain boundary magnetic field independent terms (which are expected to arize see [39, 50, 51]) have not yet 
been considered in the case of the spin chain. It is in fact not clear whether the thermodynamic Bethe Ansatz 
in its present setting is capable of providing the full boundary free energy. A less ambitious problem would be 
to consider the boundary magnetization in the framework of thermodynamic Bethe Ansatz. Indeed, as opposed 
to the full free energy, the boundary magnetization is typically given by simple integrals and not by a series of 
multiple integrals. Therefore, the comparison to our present results is probably a manageable task. We leave these 
problems for further research. 

Lastly, it would be desirable to extract the low-T behavior of the surface free energy directly out of its rep- 
resentation (3.28). For the simple and double integrals, this can be done rather straightforwardly following the 
standart method of low-T asymptotic analysis. Yet, the series of multiple integrals characterizing the function f 
(3.32) poses serious problems for extracting the low-T asymptotic analysis. We expect that all terms of the series 
will contribute to the leading orders, and thus a more refined method of low-T asymptotic analysis, possibly in the 
spirit of [28], should be implemented. 
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A Low-temperature behavior of integrals 

In this appendix, we briefly recall a lemma that allows one to carry out the asymptotic analysis of integrals 
involving the function ln[l + a(a>)] integrated versus some regular function. This lemma can be proven along the 
lines given in [35]. 

Lemma A. 1 Assume that 

• u is holomorphic in a neighborhood U ofR, 

• u has two simple zeroes 'q^' in U and %[u(A)] — > +oo with A £ U; 

2t(-l)-»+oo 
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• In 1 + e t \is holomorphic on U. 

• for any k, e ~ - TO (A~°°) uniformly in 2\ (A) — > ±00 /« {/. 

For any function f holomorphic on U with a polynomial growth along !R (A) — > ±00 in JJ, f = 0(A k ), i€N 
one /zcis 

?(+ , 2 

J /(A) In [l + e-^](U = -i J fUMA) • cU + _ ^(? ( " ) )} + 0(r 3 ) . (A.l) 



7<-> 



We stress that the O is uniform in respect to any auxiliary parameter over which / depends. If the l.h.s. in 
(A.l) is integrable in respect to this parameter so is the r.h.s. In such a case, the 0(T 3 ) remains unaffected by such 
an integration. 



References 

[1] I. Affleck, "Universal term in the free energy at a critical point and the conformal anomaly.", Phys. Rev. 
Lett. 56 (1986), 746-748. 

[2] I. Affleck, S.S. Eggert, S. Fujimoto, N. Laflorencie, and J. Sirker, "Thermodynamics of impurities in the 
anisotropic Heisenberg spin- 1/2 chain. ", J. Stat. Mech. (2008), P02015. 

[3] F.C. Alcaraz, N.M. Batchelor, R.J. Baxter, and G. R.W. Quispel, "Surface exponents of the quantum XXZ, 
Ashkin-Teller and Potts models. ", J. Phys. A: Math. Gen. 20 (1987), 6397-6409. 

[4] M. T. Batchelor and A. Klumper, "An analytic treatment of finite-size corrections in the spin-1 antiferromag- 
netic XXZ chain. ", J. Phys. A: Math; Gen. 23 (1990), L-189-195. 

[5] H. Bethe, "On the theory of metals: Eigenvalues and Eigenfunctions of a linear chain of atoms. ", Zeitschrift 
fiir Physik 71 (1931), 205-226. 

[6] F. Gohmann, M. Bortz and H. Frahm, "Surface free energy for systems with integrable boundary conditions. ", 
J. Phys. A: Math. Gen. 38 (2005), 10879-10892. 

[7] M. Bortz and J. Sirker, "Boundary susceptibility in the open XXZ-chain. ", J. Phys. A: Math. Gen. 38 (2005), 
5957. 

[8] , "The open XXZ-chain: bosonisation, Bethe Ansatz, and logarithmic corrections.", J.Stat.Mech. 

(2006), P01007. 

[9] J.L. Cardy, "Conformal invariance and universality infinite-size scaling. ", J. Phys. A: Math. Gen. 17 (1984), 
L3 85-3 87. 

[10] V.I. Cheredink, "Factorizing particles on a half-line and root systems. ", Theor. Math. Phys. 61 (1984), 977. 

[11] P. de Sa and A. M. Tsvelik, "Anisotropic spin-1/2 Heisenberg chain with open boundary conditions. ", Phys. 
Rev. B 52 (1995), 3067-3070. 

[12] P. Dorey, D. Fioravanti, C. Rim, and R. Tateo, "Integrable quantum field theory with boundaries: the exact 
g-function. ", Nucl. Phys. B696 (2004), 445^167. 



35 



[13] T.C. Dorlas, J.T. Lewis, and J.V. Pule, "The Yang-Yang thermodynamic formalism and large deviations.", 
Comm. Math. Phys. 124, 3 (1989), 365^102. 

[14] S. Eggert and S. Fujimoto, Boundary susceptibility in the spin-1/2 chain: Curie like behavior without mag- 
netic impurities, Phys. Rev. Lett. 92 (2004), 037206. 

[15] F. H.L. Essler, H. Frahm, F. Gohmann, A. Kliimper, and V. E. Korepin, "The one-dimensional Hubbard 
model. ", Cambridge University Press, 2005. 

[16] H. Frahm and A. A. Zvyagin, "The open spin chain with impurity: an exact solution. ", J. Phys.: cond. mat. 
9 (1997), no. 45, 9939. 

[17] A. Furusaki and T. Hikihara, "Boundary contributions to specific heat and susceptibility in the spin-1/2 XXZ 
chain. ", Phys. Rev. B 69 (2004), no. 9, 094429. 

[18] M. Gaudin, "Thermodynamics of a Heisenberg-Ising ring for A > 1. ", Phys. Rev. Lett 26 (1971), 1301-1304. 

[19] F. Gohmann, Private communication. 

[20] F. Gohmann, A. Kliimper, and A. Seel, "Integral representations for correlation functions of the XXZ chain 
at finite temperature. ", J. Phys. A: Math. Gen. 37 (2004), 7625-7652. 

[21] F.D.M. Haldane, "Demonstration of the "Luttinger liquid" character of Bethe-Ansatz soluble models ofl-D 
quantum fluids", Phys. Lett. A 81 (1981), 153-155. 

[22] A.G. Izergin, "Statistical sum of the six-vertex model infinite volume. ", Dokl. Akad. Nauk SSSR 297 (1987), 
331-333. 

[23] A.G. Izergin, N. Kitanine, J.M. Maillet, and V. Terras, "Spontaneous magnetization of the XXZ Heisenberg 
spin 1/2 chain. ", Nucl. Phys. B 554 (1999), 679-696. 

[24] M. Jimbo, R. Kedem, T. Kojima, H. Konno, and T. Miwa, "XXZ chain with a boundary. ", Nucl. Phys. B 441 
(1995), 437-470. 

[25] A. Kapustin and S. Skorik, "Surface excitations and surface energy of the antiferromagnetic XXZ chain by 
the Bethe Ansatz approach. ", J. Phys. A: Math. Gen 29 (1996), 1629. 

[26] N. Kitanine, K.K. Kozlowski, J.-M. Maillet, G. Niccoli, N.A. Slavnov, and V. Terras, "Correlation functions 
of the open XXZ chain I. ", J. Stat. Mech.: Th. and Exp. (2007), P10009. 

[27] N. Kitanine, K.K. Kozlowski, J.-M. Maillet, G. Niccoli, N.A. Slavnov, and V. Terras, "Correlation functions 
of the open XXZ chain II. ", J. Stat. Mech.: Th. and Exp. (2008), P07010. 

[28] N. Kitanine, K.K. Kozlowski, J.-M. Maillet, N.A. Slavnov, and V. Terras, "Algebraic Bethe Ansatz approach 
to the asymptotics behavior of correlation functions. ", J. Stat. Mech: Th. and Exp. 04 (2009), P04003. 

[29] A. Kliimper, "Thermodynamics of the anisotropic spin-1/2 Heisenberg chain and related quantum chains. ", 
Zeit. Fur Phys.: Cond. Mat. B 91 (1993), 507-519. 

[30] A. Kliimper, M. T. Batchelor, and P. A. Pearce, "Central charges for the 6- and 19-vertex models with twisted 
boundary conditions. ", J. Phys. A: Math. Gen. 24 (1991), 3111-3133. 



36 



[31] A. Kliimper, T. Wehner, and J. Zittartz, "Conformal spectrum of the 6-vertex model. ", J. Phys. A: Math. Gen. 
26(1993), 2815-2827. 

[32] T. Koma, "Thermal Bethe-Ansatz method for the one-dimensional Heisenberg model. ", Prog. Theor. Phys. 
78(1987), 1213-1218. 

[33] , "Thermal Bethe-Ansatz method for the spin-Xjl XXZ Heisenberg model", Prog. Theor. Phys. 81 

(1989), 783-809. 

[34] V. E. Korepin, "Calculation of norms ofBethe wave-functions. ", Comm. Math.Phys. 86 (1982), 391. 

[35] K. K. Kozlowski, "Low-T asymptotic expansion of the solution to the Yang-Yang equation. ", math-ph: 
11126199. 

[36] K. K. Kozlowski, J.-M. Maillet, and N. A. Slavnov, "Low-temperature limit of the long-distance asymptotics 
in the non-linear Schrodinger model. ", J.Stat.Mech. (2011), P03019. 

[37] E.H. Lieb and W. Liniger, "Exact analysis of an interacting Bose gas. I. The general solution and the ground 
state. ", Phys. Rev. 130 (1963), 1605-1616. 

[38] B. M. McCoy, "Spin correlation functions in the XY model. ", Phys. Rev. 173 (1968), 531-541. 

[39] B. Pozsgay, On 0(1) contributions to the free energy in Bethe Ansatz systems: the exact g-function, JHEP 08 
(2010), 090. 

[40] D. Ruelle, "Statistical mechanics: rigorous results", W.A. Benjamin, Inc., 1969. 

[41] E. K. Sklyanin, "Boundary conditions for integrable quantum systems.", J. Phys. A: Math. Gen. 28 (1988), 
2375-2389. 

[42] J. Suzuki, "Relationship between d-dimensional quantal spin systems and (d+\)-dimensional Ising systems. ", 
Prog. Theor. Phys. 56 (1976), 1454-1469. 

[43] M. Suzuki, "Transfer-matrix methods and Monte Carlo simulation in quantum spin systems. ", Phys. Rev. B 
31 (1985), 2957-2965. 

[44] M. Takahashi, "Thermodynamics of one-dimensional solvable models. ", Cambridge University Press, 1999. 

[45] M. Takahashi, "One-dimensional Heisenberg model at finite temperature.", Prog. Theor. Phys. 42 (1971), 
1289. 

[46] , "Correlation length and free energy of the S= 1/2 XXZ chain in a magnetic field. ", Phys. Rev. B 44 

(1991), 12382-12394. 

[47] , "Correlation length and free energy of the S= 1/2 XYZ chain. ", Phys. Rev. B 43 (1991), 5788-5797. 

[48] C. A. Tracy and H. Vaidya, "One particle reduced density matrix of impenetrable bosons in one dimension 
at zero temperature. ", J. Math. Phys. 20 (1979), 2291-2312. 

[49] O. Tsuchiya, "Determinant formula for the six-vertex model with reflecting end. ", J. Phys. A: Math. Gen. 39 
(1998), 5946-5951. 



37 



[50] F. Woynarovich, "0(1) contribution of saddle-point fluctuations to the free energy of Bethe Ansatz systems. ", 
Nucl. Phys. B700 (2004), 331. 

[51] , "On the normalization of the partition function of Bethe Ansatz systems. ", Nuclear Physics B 852 

(2011), 269-286. 

[52] T. T. Wu, "Theory of Toeplitz determinants and the spin correlation functions of the two-dimensional Ising 
model I. ", Phys. Rev. 149 (1966), 380-401. 

[53] C. N. Yang and C. P. Yang, "Thermodynamics of a one-dimensional system of bosons with repulsive delta- 
interactions. ", J. Math. Phys. 10 (1969), 1115-1122. 

[54] A. A. Zvyagin and A. V. Makarova, "Bethe-Ansatz study of the low -temperature thermodynamics of an open 
Heisenberg chain", Phys. Rev. B 69 (2004), no. 21, 214430. 



38 



